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Abstract. Ultracold neutral bosons in a rapidly rotating atomic trap have been 
predicted to exhibit fractional quantum Hall-like states. We describe how the 
composite fermion theory, used in the description of the fractional quantum Hall effect 
for electrons, can be applied to interacting bosons. Numerical evidence supporting the 
formation of composite fermions, each being the bound state of a boson and one flux 
quantum, is shown for filling fractions of the type v = p/(p + 1), both by spectral 
analysis and by direct comparison with trial wave functions. The rapidly rotating 
system of two-dimensional bosons thus constitutes an interesting example of "statistical 
transmutation," with bosons behaving like composite fermions. We also describe the 
difference between the electronic and the bosonic cases when p approaches infinity. 
Residual interactions between composite fermions are attractive in this limit, resulting 
in a paired composite-fermion state described by the Moore-Read wave function. 
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1. Introduction 

Under appropriate circumstances, strong short range repulsion between bosons can 
mimic the Pauli exclusion principle of fermions. This has been known for a long time in 
the case of bosons in one dimension pQ, and the Lieb-Liniger criterion |2] discriminates in 
which case this fermionization occurs. Experiments based on Bose-Einstein condensates 
of atomic gases allow a test of the concept of statistical transmutation of interacting 
bosons into fermions. Such behavior has been recently observed in experiments j31 E] ■ 

Here we discuss mapping of bosons into a different kind of fermions, called composite 
fermions. Each composite fermion is the bound state of a boson and one quantized 
vortex of the wave function, often conveniently though of as the bound state of a boson 
and one flux quantum, O — hc/e. Composite fermions are expected to occur when 
ultracold bosons are placed in a rapidly rotating atomic trap. Such rotation leads to 
a variety of rich phenomena, such as the Abrikosov lattice of vortices jSl ©• For high 
rotation frequency, it has been predicted [3 El El HH] that bosons, under appropriate 
conditions, will behave like two dimensional electron systems (2DES) in the fractional 
quantum Hall effect (FQHE) regime. [More precisely: A strong confinement along 
the rotation axis is required so the system can be considered as two dimensional. A 
fine tuning between the rotation frequency and the harmonic trap frequency must be 
reached so that the centrifugal force compensate the harmonic trapping force. Further, 
the interaction between bosons has to be weak enough so that any mixing with higher- 
lying Landau levels is negligible. Thus the system of neutral bosons subjected to a 
Coriolis force is equivalent to a system of charged particles in an external magnetic 
field. The filling factor v is the ratio of the number N of particles over the number Ny 
of quantized vortices that would be present in the system if it were a Bose condensate.] 
This system of bosons in the lowest Landau level (LLL) is an ideal playground where 
theories of the FQHE can be tested, and extended into new regimes. In some respects, 
the bosonic system offers more flexibility than the electronic one. The strength of the 
interaction can be tuned using the Feschbach resonance, and there are many choices for 
the condensed atoms; that allows the possibility of studying systems with different types 
of interactions (delta function [7j, hollow core [TT] or dipole-dipole [T^ ) and quantum 
statistics (Bose or Fermi). 

For the electron system in the lowest Landau level, Laughlin's wave functionfTBl 
provides a good approximation for the ground state at filling factor 1/3. However, a 
larger class of fractions of the type pj (2p± 1) appear in experiments, and are understood 
through the formation of composite fermions (CFs) [H], namely electrons bound to an 
even number of quantized vortices. The problem of interacting electrons maps onto 
free composite fermions in an effective magnetic field, and the fractions p/(2p ± 1) are 
explained as the integral quantum Hall effect of composite fermions. 

In this paper we focus on the case of bosonic atoms with only one hyperfine species 
(spinless bosons). In this case, the interaction between atoms occurs through s-wave 
scattering. Such a system has the property that the Laughlin wave function |13j is 
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the exact ground state at filling factor v = 1/2 7J. However, one can ask: Do 
composite fermions also occur in rotating two-dimensional bose systems? The composite 
fermionization of the bosonic system will have many testable consequences beyond 
Laughlin's wave function, most important being the prediction of incompressible states 
at many fractions of the type v = p/(p ± 1). Such fractions can also be deduced 
from the fermion-boson mapping based on a gauge field approach [13 HEl EI- We 
show here numerical evidence of the accuracy of the CF description by identifying 
the incompressible states, by spectral analysis, and by direct comparison of the Jain's 
CF wave functions with exact wave functions. We shall use the spherical geometry 
to establish detailed and quantitative correspondence between interacting bosons and 
free fermions. This study extends previous works on disk |18| |SJ E2] and toroidal 
geometries jHj- An advantage of the spherical geometry over the disk geometry is in 
facilitating the identification of incompressible states of the type 2/3 and 3/4. 

In section |2l we present the effective Hamiltonian that describes the system. Section 
0] is devoted to the Laughlin wave function, which is the exact ground state of the 
rotating bosons interacting through s-wave scattering at filling factor v = 1/2. We 
present the spherical geometry in section which provides a convenient framework 
for numerical calculations. In section we give an overview of the CF theory and 
its extension to bosons. Section |H1 shows numerical evidence for the presence of 
incompressible states at several fillings, and demonstrates the suitability of the CF 
theory for their description. Finally in section [JJ we discuss the state at v — 1, which 
requires a consideration of pairing of composite fermions, described by the Moore- Read's 
Pfaffian wave function. 

2. Effective Hamiltonian 

In the rotating frame, the Hamiltonian that describes N atoms of mass m in a rotating 
harmonic trap is given by : 

( 2 2\ / 2 i 2\ i 2 2 

\u) r -u) ) [x\ +y i ) + u z z i 
+ *£V(r i -r j ), (1) 

i<j 

where u r is the radial harmonic trap frequency, uj z is the axial harmonic trap frequency 
and u the trap rotation frequency (rotation along the z axis). For ultracold bosonic 
gases, the interaction is dominated by s-wave scattering. Thus the interaction can be 
approximated by a delta function : 

V (r) = ^^l 5 m ( r ) , (2) 
m 

with a s the s-wave scattering length. 

We will now assume the critical rotation frequency u = u r . Thus in the plane 
perpendicular to the rotation axis, the system is equivalent to charged particles in a 
magnetic field along the z axis with cyclotron frequency u c = 2u. We also require 



H = L (Pi - x r<) + - 
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strong confinement along the z axis, so the wave function in this direction is the ground 
state of the harmonic oscillator along z. The effective 2D interaction is given by : 



V (r) = gl 2 5 {2) (r) with g = V32^ huy-, (3) 



where r is now the particle coordinate in the x — y plane, l z = Jh/mui z is the 



characteristic length of the z axis oscillator and I = yfi/2muj is the magnetic length. 
The effective Hamiltonian can be then written as : 

H 2d = £ 7T (Pi + 9l 2 E J(2) ( r * - r i) • ( 4 ) 

8=1 lm i<j 

This two dimensional regime is formally equivalent to the FQHE for 2DES. The filling 
factor v is defined as v = hp/2muj where p is the areal boson density. It can be also 
expressed as the ratio N/Ny with Ny the number of quantized vortices which is the 
analog of the number of flux quanta in 2DES systems. We will assume in the 
following that the interaction strength is small enough so that Landau level mixing can 
be neglected, and all bosons are in the LLL. The effective Hamiltonian is then reduced 
to the interaction term alone : 

tfLLL^E^" 1 ;)- (5) 

i<j 

Note that due to the bosonic statistics, we can have filling factors greater that one and 
still stay entirely in the LLL. Such a regime is forbidden in the fermionic case. 

3. v = 1/2: the Laughlin state 

In the LLL, the one-body eigenfunctions in the symmetric gauge are given by : 

Mz)= m * , fflVl'W (6) 



V27r/ 2 2 m m! W. 

with z = x + iy, the positive integer m is the angular momentum of the state (L z = km) 
and / is the magnetic length. Any N-body wave function for particles in the LLL can 
be expressed as a polynomial V of the Z{ coordinates up to a global Gaussian factor : 

*( Zl ,...,z N ) =V(z u ...,z N ) e-^NV*". (7) 

The Laughlin wave function ^3] was introduced in the case of 2DES to explain the 
fractional quantum Hall state at filling factor v = 1/3. It describes a droplet of 
incompressible fluid with a mean density corresponding to this fraction. It can also 
be generalized for bosons : 

laughlin (ZX,...,Z N ) = n>-**) 2 e-^N 2 /^. (8) 

In the case of 2DES, the Laughlin wave function is a very good approximation for the 
true ground state. But for the system we are considering, this is the exact zero energy 
ground state of the effective Hamiltonian Eq.(|5|) at filling factor v = 1/2. It is also 
the zero energy ground state with the smallest total angular momentum (equal to the 
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highest filling factor). An arbitrary function built from the product of ^Laughlin an d 
any symmetric polynomial of the Zi coordinates is also a zero energy ground state but 
has higher angular momentum and thus corresponds to a filling factor v < 1/2. These 
states represent gapless edge excitations or quasihole excitations of the 1/2 system. 
Quasielectrons are nucleated from the Laughlin state when reducing the number of flux 
quanta (or the angular momentum). Those have a nonzero gap in the thermodynamical 
limit JU|- These properties can be seen in exact diagonalizations. Figure ^ displays 
the spectrum for seven particles. The first zero energy eigenstate appears for angular 
momentum L z =42 which is in agreement with Eq. |H] (L z = N(N — 1)). Higher angular 
momentum values also exhibits zero energy eigenstates as expected from the above 
discussion. 
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Figure 1. Lower part of the energy spectrum for 7 bosons on the disk geometry. The 
Laughlin state corresponds to the zero energy ground state at L Z =A2. 



4. The spherical geometry 

The planar, or disk, geometry is natural for the system we are dealing with. 
Unfortunately this geometry is plagued by edge effects when doing numerical 
calculations. Because we are interested in the bulk properties of our system, it is more 
suitable to work with compact geometries, such as a torus or a sphere. In this paper 
we will use the last one, which is easier to handle due to its SU(2) symmetry, enabling 
a classification of the states with respect to their total angular momentum (denoted by 
L). Let us introduce the spinor coordinates : 

Uj = cos (0j/2) e^ /2 and v, = sin (9j/2) e"^ /2 . (9) 
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The spherical geometry requires us to introduce a magnetic monopole at the center of 
the sphere. Its radius R is then related to the number of flux quanta that pierce it : 

R = l^j Ny/2. (10) 

Solutions of the one-body problem are given by the monopole harmonics [2D| (a 
generalization of the spherical harmonics), which take the following form in the LLL : 



(25 + 1)! . +m ^ 



« i+ V" m , 11 



\|47r(S-m)!(S + m) 

where m is the projection of the angular momentum along the z axis. The angular 
momentum quantum number S for the LLL is related to Ny by the relation 25 = Ny. 
In a given Landau level, the two-body interaction can be parameterized using a set 
of 25 + 1 number {V m } called the pseudo-potentials where m is related to the 
relative angular momentum between the two particles. For spinless bosons, only even- 
m potentials are relevant. In this formalism, the delta function interaction corresponds 
to all V m = except Vq. Longer range interaction can be simulated by turning on 
the other pseudo-potentials, the first relevant one for bosons being Vi. For trial wave 
functions such as the Laughlin wave function, we can pass from the disk geometry to 
the spherical geometry using the stereographic projection (221 • F rom a practical point 
of view, we drop the Gaussian factor and make the replacement : 



3) 



(uiVj-UjVi). (12) 



Jain's wave functions are conveniently constructed directly in the spherical geometry. 



5. Composite fermions 

In the case of 2DES, the CF theory ^3] maps the problem of strongly interacting 
electrons in the LLL onto weakly interacting fermions called composite fermions. Each 
of these fermions behaves, in a sense, as if it has captured an even number q of flux 
quanta. Thus, the CFs feel a reduced number of external flux quanta Nl = N<p — qN, 
and the effective filling factor v* for this system is related to the true filling factor v by 
the relation : 

z/ = — ^ . (13) 

1 + qv* v ' 

When the CFs completely occupy an integer number of pseudo-Landau levels, i.e., when 
v* is an integer p, they produce incompressible states at the Jain's principal sequences 
v = p/(qp + 1), which includes the Laughlin states at v — 1/(5+1). The CF theory 
also guides construction of trial wave functions. They are of the form = 
where <t> v * is the Slater determinant corresponding to the uncoupled CFs and 4>i is the 
filled Landau level (i.e. the Jastrow factor) : 

0i = n & - . ( i4 ) 

i<j 
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which encodes the flux attachment. Because the product (fil4> u * has components in higher 
Landau levels, it is projected onto the LLL, as indicated by the projection operator Plll- 
The CF construction trivially generalizes for bosons, with the slight change that 
q is now taken as an odd integer. The CFs built from this process obey fermionic 
statistics which means we can use the same construction as in the electronic case 
(i.e. filling pseudo Landau levels). So for bosons, the principal sequence occurs at : 

^ = ^r, (15) 

and the corresponding trial wave functions are = Trial wave functions 

can also be written for the spherical geometry j2H • For the ground state associated 
to the fraction v — p/(p + 1), the relation between S and iV is given by 

2S =^^N -p-1. (16) 
p 

In general, the trial wave functions are more complicated than Laughlin's, but many 
of their essential properties can be explicitly derived, such as the maximum angular 
momentum of the excitons, or that the ground state has L = total angular momentum 
(which follows from the fact that the ground state contains filled shells of composite 
fermions) . 



6. Evidence from exact diagonalization 

6.1. spectrum analysis 

Exact diagonalization has been widely used to investigate the FQHE since its early 
days. Even for a small number of particles, clear signatures of FQHE are exhibited in 
numerical data [26J. In our calculations, we use the spherical geometry and diagonalize 
the Hamiltonian in a fixed L z subspace. Because we are mainly interested in the low 
energy physics, numerical diagonalization is performed using the Lanczos algorithm, 
which allows a determination of eigenvectors and eigenvalues for large matrices. The 
largest Hilbert space dimensions we can handle are of the order of 10 7 , which enables 
us to deal with up to iV ~ 15 — 30 particles, depending on the filling factor. Candidates 
for incompressible states are identifiable by their spectrum, because they have a clearly 
non-degenerate ground state (with the total orbital angular momentum L — 0). FigureEl 
shows the spectrum of a typical such spectrum. Once all reachable spectra have been 
evaluated, we start searching for sets of (25*, N) values that satisfy a linear relation of 
the type : 

2S = —N — shift. (17) 
v 

If all (2S, N) obeying this relation are incompressible, and the gap extrapolates to 
a non-zero value in the thermodynamic limit, then the corresponding fraction should 
exhibit FQHE in the real system. We find this to be the case for fractions v = 1/2, 
2/3 and 3/4 JU]- For these three fractions, the relation between 2S and N is the one 
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predicted by Eg. ([16]) (including the non-trivial shift). This is a strong hint that the 
CF theory correctly describes these fractions. Other spectral evidences such as the 
maximum angular momentum of collective mode, or the L value of the guasiparticule 
excitations, further reinforce the CF hypothesis. But the most direct and compelling 
verification is provided by the overlap calculations. 




2 4 6 8 10 12 14 5 10 15 20 

L 2S 



Figure 2. Left panel : The low energy spectrum of a typical incompressible state in 
the spherical geometry (here for N = 9 at v = 1/2). The ground state at L = (in 
this case the Laughlin state) is clearly separated from the excitations. Right panel : 
The candidates for incompressible states are marked on the N — 2S plane, where N 
is the number of bosons and 2S the number of flux quanta through the surface of the 
sphere. Each fraction corresponds to a line with a non-trivial shift (plotted here for 
fractions v — 1/2, v — 2/3 and v = 3/4). 



6.2. Overlap calculations 

The strength of Jain's CF theory is that it predicts not only the filling factors for 
incompressibility but also provides explicit wave functions for ground and excited states. 
Numerical comparisons can be carried out between these wave functions and those 
obtained from exact diagonalization. The overlap between the exact state ty ex and the 
trial wave function \l> trial is defined as follows : 



O = 





Atrial) | 2 


1 


I (Atrial |^trial)| 



(18) 



Scalar products are calculated using Monte Carlo integration. This requires an 
evaluation of the exact state in the real space, involving a computation of a huge 
number of permanents, which is a time consuming task (see j2H] for more details). 
Results of overlap calculations [23] are displayed in tabled for fractions v = 1/2 and 
2/3 and for both ground and excited states. They prove that the CF theory is in good 
agreement with exact diagonalizations. The CF description is also robust with respect 
to changes of the interaction. Figure El shows the evolution of the ground state when we 
tune the interaction from pure delta function to a longer range interaction by adding 
some Vi pseudo-potential term. Obviously the CF approach is valid for a wide range of 
interactions (including the Coulomb interaction j2Sj)- For other fractions of the principal 
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Table 1. Overlaps between exact wave functions and trial wave functions for ground 
and excited states for various system sizes, denoted by O gr and O ex . For v = 1/2, 
the trial wave function for the ground state is Laughlin's wave function, and the trial 
wave functions for the first excited states at angular momentum L are Jain's wave 
functions. For v — 2/3, both the ground state and the first excited state overlaps 
involve the Jain's wave functions. The statistical uncertainty in the last two digits 
from the Monte Carlo integration is shown in parentheses when it is larger than 10 -5 . 



v N 


S r 


L 




v N 




L 




1/2 4 


1.0000 


4 


0.9972 


2/3 4 


1.0000 


2 


1.0000 


5 


1.0000 


4 


0.9965 


6 


0.9850 


4 


0.7544(05) 


6 


1.0000 


■5 


0.9959 


8 


0.9820(10) 


5 


0.8701(14) 


7 


1.0000 


5 


0.9954 


10 


0.9724(89) 


G 


0.855(12) 


8 


1.0000 


6 


0.9945 










9 


1.0000 


6 


0.9954(2) 











sequence p/(p + 1), the situation becomes worse as p — > oo, i.e. as v approaches 1. A 
more careful study of this limit remains to be done. 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 



v 2 /v„ v 2 /v 

Figure 3. Overlap of the exact ground state of an effective Hamiltonian, with an 
additional longer range component V2 in addition to the contact interaction Vo, with 
Laughlin's and Jain's trial wave functions at v — 1/2 (left panel) and v = 2/3 (right 
panel), as a function of V2/V0. The overlaps in the left panel are evaluated exactly, 
and those in the right panel by Monte Carlo method (statistical error not shown). 

7. v — 1: the PfafRan state 

The 2DES does not show FQHE at v = 1/2 |2Zj. It is the accumulation point of 
the Jain's principal sequence v = p/(2p + 1), and exhibits behavior akin to a Fermi 
sea. This property can be understood naturally within the CF theory, as the effective 
magnetic field vanishes at v = 1/2. Thus CFs behave like free fermions in a zeroth-order 
approximation. 

The obvious question is then whether a Fermi sea of composite fermions occurs 
for bosons at the corresponding filling factor v — 1. This would be true if the model 
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of free composite fermions is qualitatively valid here. It has been argued the residual 
interaction can lead to pairing of CFs[2H]- This proposal has been introduced to explain 
the incompressible state observed at v = 5/2 in 2DES. Moore and Read[2H] have 
proposed an explicit trial wave function to describe the paired ground state which is 
analogous to the BCS wave function in real space for a given number of particles. The 
corresponding wave function for bosons can be written for the spherical geometry as : 

1 



*Pf (zi, 



Pf 



UiVj-UjVi, iKj 



I I (UiVj - UjVi) , 



(19) 



where Pf(A) is the Pfaffian of the skew-symmetric N x N matrix A (N even), defined 
as : 

Pi (A) = 2Ze CT A CT (i) -( 2 )vi (T (3) .(4)...yi (T (jv-i)<T(Ar), (20) 



where the sum runs over all permutations and e a is the signature of the permutation. 
If such a state is present, it should appear for even numbers of particles at 2S = N — 2. 
Exact diagonalization clearly shows incompressible states at these values of (2S, N) and 
also parity effect (see figure EJ). Thermodynamic extrapolation of the gap leads to a 
non-zero value [TUj . 
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(a) N = 12 2S = 10 
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(c) N = 14 2S = 12 
12 



14 



Figure 4. Spectra at rilling factor v = 1 for N = 12,25 = 10 (a), N = 13,25 = 11 
(b) and N = 14, 25 = 12 (c) in the spherical geometry. The incompressible states only 
appear for even numbers of particles. The relation between the number of bosons and 
flux quanta for these incompressible states is the one predicted by the Moore-Read 
wave function. 

Direct comparison between the exact ground state and the Moore- Read (or Pfaffian) 
wave function can also be performed j2S|- The latter can be generated by exact 
diagonalization using the fact that Eq. (fT5j) i s the exact zero-energy ground state of 
the three-body hardcore Hamiltonian [HOI IS] '■ 

Hpf = ^ (2) ( u i v i ~ u J v i) ^ (2) ( u i v k - U kVj) ■ (21) 

i<j<k 

The overlaps are thus obtained by a simple scalar product of two ground state vectors. 
Overlap results in table El show that the description based on the Pfaffian wave function 
is valid for the delta function interaction; figure El demonstrates its robustness for longer 
range interactions as well. 
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Table 2. Overlap between the exact ground state at 2S = N — 2 and the Pfaffian 
state for different system sizes. 



N 


4 


6 


8 


10 


12 


14 


16 


O 


1.0 


0.9728 


0.9669 


0.9592 


0.8844 


0.8858 


0.8833 



1 

0.8 

1) 0.6 
> 
a 

0.4 
0.2 

0.1 0.2 0.3 0.4 0.5 

v 2 /v 

Figure 5. Overlap of the Pfaffian state with the exact ground state of the effective 
Hamiltonian with an additional longer range component Vi- The overlaps are 
calculated for N = 14 and N — 12 bosons. The description in terms of the Pfaffian 
state is robust with respect to long range interaction up to a critical value of V2/V0. 

8. Conclusion 

Both analytical and numerical results suggest that ultracold boson atoms in rotating trap 
should display a very rich physics in the low filling factor regime. The analogy with the 
FQHE in 2DES suggests the use of the CF approach for the bosonic system, and indeed, 
numerical tests prove its validity. The robustness of the CF description for fractions 
v = 1/2 and v = 2/3 is clearly demonstrated. Nevertheless, there are some major 
differences from the 2DES systems. The residual interactions between the composite 
fermions tend to complicate the description of fractions p/(p + 1) with increasing p, 
leading to a paired state at v — 1 well described by a BCS-like Pfaffian wave function. 
Also, in this paper, we do not address the v > 1 case. CF theory allows construction 
of the sequence v = p/{p — 1) by reversing the attached flux quanta. But the overlap 
calculations tend to indicate that this description is not valid (at least for v = 2 and 
v = 3/2). Other states have been predicted |Hj to appear at filling factors v = k/2 (k 
being an integer value), correspond to Read-Rezayi states. They are generalizations of 
the Pfaffian state and also exhibits excitations obeying non-Abelian statistics. But it is 
not yet clear whether this description is valid in the thermodynamic limit for the delta 
function interaction ^Uj- Numerical evidence [22] indicates that longer range interaction 
(like dipole-dipole interaction) help to stabilize such exotic states. 
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